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Refinement  of  a  Semi-Empirical  Model  for  the  Microwave  Emissivity  of 
the  Sea  Surface  as  a  Function  of  Wind  Speed.  (December  1995) 

David  Jacob  Kohn,  B.S.,  Michigan  Technological  University 
Chair  of  Advisory  Committee;  Dr.  Thomas  T.  Wilheit 

In  1979,  Wilheit  introduced  a  sea  surface  emissivity  model  for  microwave 
frequencies.  This  model  is  used  in  a  radiative  transfer  model  (RTM)  to  obtain  simulated 
brightness  temperatures  for  various  atmospheric  conditions.  The  brighmess  temperatures 
are  used  to  obtain  algorithms  for  atmospheric  variables.  These  algorithms  can  then  be  used 
to  retrieve  atmospheric  variables  from  the  microwave  measurements.  Therefore,  it  is 
important  to  get  the  surface  emissivity  right. 

Several  changes  are  made  to  Wilheit's  sea  surface  emissivity  model.  The  first 
change  to  the  model  is  to  the  model's  treatment  of  multiple  reflections.  Multiple  reflections 
are  now  treated  as  if  the  radiation  is  reflected  back  into  the  view  path  of  the  microwave 
sensor.  This  change  lowered  the  computed  emissivity  of  the  sea  surface;  which  is  more 
representative  of  observations  without  sea  foam.  The  second  change  is  made  to  the  sea 
surface  roughness  parameter.  An  increase  in  roughness  is  needed  at  frequencies  above  16.6 
GHz  and  a  decrease  below  16.6  GHz.  The  roughness  is  increased  to  132%  of  the  Cox  and 
Munk  roughness  at  37  GHz  and  30%  of  the  roughness  at  6.6  GHz.  The  last  change  to  the 
model  is  in  the  treatment  of  sea  foam.  The  foam  effect  is  now  a  smooth  transition  of 
increasing  foam  as  the  wind  speed  increases;  instead  of  being  switched  on  at  7  ra/s.  These 
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changes  yield  an  improved  sea  surface  model  that  is  within  2%  of  the  sea  surface  emissivity 
given  by  the  Wentz  sea  surface  emissivity  functions. 
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1.  INTRODUCTION 

Surface  observations  are  needed  to  initialize  forecast  models.  As  forecast  models 
become  more  sophisticated,  accuracy  and  increased  spatial  resolution  of  observations 
become  more  important.  Since  the  majority  of  the  earth's  surface  is  ocean,  a  large  part  of 
the  earth  is  not  covered  by  surface  observations.  These  observations  are  needed  for  an 
accurate  initialization  field  in  forecast  models.  Satellite  remote  sensing  can  provide  the 
necessary  atmospheric  parameters  over  the  ocean  to  fill  this  gap  in  the  surface  observation 
network. 

The  Multi-frequency  Imaging  Microwave  Radiometer  (MIMR),  a  passive 
microwave  sensor,  is  planned  to  be  launched  in  the  year  2000.  One  will  be  able  to  retrieve 
general  atmospheric  parameters  from  the  MIMR  frequencies.  A  radiative  transfer  model 
(RTM)  will  be  used  to  simulate  the  expected  brightness  temperatures  in  various 
atmospheres  for  the  MIMR  frequencies.  These  simulated  brightness  temperatures  will  be 
regressed  against  the  various  atmospheric  parameters  to  obtain  operational  algorithms. 

Currently,  the  Special  Sensor  Microwave/Imager  (SSM/I)  is  a  remote  sensing 
instrument  used  for  data  collection.  The  SSM/I  is  flown  on  polar  orbiting  satellites  and  has 
global  coverage  approximately  twice  daily.  This  sensor  measures  upwelling  radiation  from 
the  earth's  surface.  However,  a  uniform  reflective  surface  is  needed  to  make  accurate 
inferences  of  atmospheric  variables.  Therefore,  for  many  parameters  the  only  useful  data 
are  collected  over  oceanic  areas. 


The  journal  model  for  this  thesis  is  UtiQ  Journal  of  Applied  Meteorology. 
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The  SSM/I  radiances  are  reported  as  brightness  temperatures.  Brightness 
temperatures  are  measured  in  the  horizontal  and  vertical  polarizations  at  19.35,  37.0,  and 
85.5  GHz.  The  22.235  GHz  channel  measures  only  the  vertical  polarization.  Brightness 
temperatures  are  a  combination  of  surface  emission,  upward  radiation  by  atmospheric 
constituents,  and  downward  emission  and  cosmic  background  radiation  reflected  off  the  sea 
surface  (Fig.  1).  Various  geophysical  parameters;  cloud  liquid  water  (CLW),  precipitable 
water  (PW),  and  sea  surface  wind  speed  (WS);  can  be  retrieved  from  combinations  of 
brightness  temperatures.  Also,  sea  surface  temperatures  (SST)  are  needed  to  initialize  the 
forecast  models  which  cannot  be  obtained  from  the  SSM/I  frequencies.  However,  the 
MIMR  having  additional  frequencies  will  provide  all  the  above  parameters. 

An  integral  part  of  the  RTM  is  the  sea  surface  emissivity  calculation.  Wind 
roughens  the  sea  surface  and  produces  foam.  This  roughness  and  foam  cause  changes  in 
sea  surface  emissivity  and  therefore  in  the  radiances.  Wentz  has  developed  two  empirical 
sea  surface  emissivity  functions.  He  developed  one  for  the  scanning  multichannel 
microwave  radiometer  (SMMR)  frequencies  and  the  other  for  the  SSM/I  frequencies. 
Although,  these  functions  work  well  at  describing  the  surface  emissivity,  they  are  limited  to 
the  frequencies  and  incidence  angles  for  which  they  were  developed.  A  general  model  for 
the  surface  emissivities  is  needed  to  be  able  to  simulate  brightness  temperatures  at  other 
frequencies  and  incidence  angles.  The  objective  of  this  research  is  to  refine  the  surface 
emissivity  model  to  better  reflect  the  relationship  between  wind  speed  and  sea  surface 


emissivity. 


Fig.  1.  Two  component  (single  layer  atmosphere  and  surface)  radiative  transfer 
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2.  DESCRIPTION  OF  THE  SEA  SURFACE  EMISSIVITY 
FUNCTIONS  AND  MODEL 

a.  Scanning  multichannel  microwave  radiometer  (SMMR) 

From  the  scanning  multichannel  microwave  radiometer  (SMMR)  sensor  data, 
Wentz  (1983)  developed  an  expression  to  retrieve  sea  surface  emissivity.  He  separated  the 
sea  surface  emissivity  into  two  parts.  The  first  part  is  the  contribution  from  a  smooth  sea 
surface.  The  second  part  is  the  contribution  from  wind  induced  roughness  and  foam.  The 
total  emissivity  (E)  can  be  written  as 

E=Es+Er  (1) 

where  is  the  smooth  surface  emissivity  and  E^  is  the  rough  surface  and  foam  part  of  the 
emissivity. 

The  smooth  surface  contribution  to  emissivity  is  computed  from  the  dielectric 
constant  of  sea  water  in  conjunction  with  the  Fresnel  relations  (Smith  1988).  Klein  and 
Swift  (1977)  reviewed  experimental  data,  in  the  1  to  4  GHz  range,  and  derived  an 
expression  for  the  dielectric  constant  of  sea  water  as  a  function  of  temperature,  salinity, 
conductivity,  frequency,  and  polarization.  Wentz  (1983)  used  this  function  for  the 
dielectric  constant  and  the  Fresnel  relation  to  derive  an  expression  for  E^.  The  resulting 
expression  follows 


Es  = 


So+Si  T+S2  +54(0,-49'') 

Ts 


(2) 


where 


r=7’,-273.16 


(3) 
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and  T,  is  the  surface  temperature  in  Kelvin  (K).  This  expression  only  applies  to  incidence 
angles  of  49®  ±0.5".  The  coefficients  for  the  expression  at  SMMR  frequencies  are  in 
Table  1.  These  coefficients  were  found  by  using  a  least  squares  fit  to  the  available  data 

(Wentz  1983). 

After  subtracting  the  brightness  temperature  (TJ  attributed  to  the  smooth  surface 
effect  from  the  SMMR  data,  the  remainder  is  considered  to  be  due  to  a  wind  roughened  sea 
surface  and  sea  foam.  The  amount  of  roughening,  as  the  wind  blows  across  the  sea  surface, 
is  related  to  the  friction  velocity  (u.)  (Ross  and  Cardone  1974).  In  an  unstable  boundary 
layer,  cold  air  over  warm  water,  there  is  more  wind  stress  (u.)  for  a  given  wind  speed  at  a 
height  over  the  water  than  in  a  neutral  or  stable  boundary  layer.  The  remaining  brighmess 
temperature  is  regressed  with  u.  to  derive  an  expression  for  the  rough  surface  and  foam 
contribution  to  the  emissivity  as  a  function  of  u.  (Wentz  1983).  Wentz  et  al.  (1986)  revised 
the  u»  relation  to  one  that  uses  wind  speed  in  m/s  at  19.5  meters.  The  new  expressions  for 
the  rough  surface  and  foam  contribution  (E^)  to  the  emissivity  are 


Er  =  m\WS 

WS  <  Imis 

(4) 

Er  =  mxWS+Q.OS{jni-mi){WS-l? 

7<WS<nm/s 

(5) 

Er  =  ni2WS—  12(/W2  —mi) 

WS>  \lmls 

(6) 

where  WS  is  the  wind  speed.  The  m  coefficients  are  in  Table  2.  The  coefficients  were 
found  by  a  least  squares  fit  to  the  residuals  of  the  brightness  temperature  (Wentz  et  al. 


1986). 
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TABLE  1.  Regression  coefficients  for  smooth  surface  brightness  temperature 
(Wentz  1983). 


SMMR 

So 

s. 

^2 

Ss 

S4 

Channel 

(K) 

(K-') 

(K-^) 

(K/deg) 

6.6  V 

1.3759  X  10^ 

2.368  X  10-' 

1.565  X  10-" 

-2.311  X  10-" 

2.03 

6.6  H 

0.7107  X  10" 

0.891  X  10' 

1.000  X  10" 

-2.497  X  10"* 

-1.28 

10.7  V 

1.4452  X  10" 

-0.336  X  10  ' 

2.076  X  10" 

-2.497  X  10^ 

2.05 

10.7  H 

0.7559  X  10" 

-0.935  X  10-' 

1.371  X  10" 

-1.661  X  lO"* 

-1.32 

18  V 

1.5750  X  10" 

-3.936x  10' 

-2.048  X  10-^ 

2.08 

18  H 

0.8444  X  10" 

-3.675x  10' 

-1.40 

21V 

1.6252x10" 

-4.916x10’ 

-1.775  X  lO"" 

2.10 

21  H 

0.8802  X  10" 

-4.546x  10' 

1.699  X  10" 

-1.477  X  lO"* 

-1.43 

37  V 

1.8493  X  10" 

-7.405  X  10  ' 

1.694  X  10" 

-0.539  X  lO"* 

2.11 

37  H 

1.0524  X  10" 

-7.666  X  10  ' 

1.718  X  10  " 

-1.033  X  10-" 

-1.59 
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b.  Special  sensor  microwave  imager  (SSM/I) 

Wentz  (1992)  also  developed  an  expression  for  sea  surface  emissivity  using  SSM/I 
frequencies  which  is  similar  to  his  SMMR  sea  surface  emissivity  expression.  Wentz' 
SSM/I  sea  surface  emissivity  function  can  be  summarized  as: 

E=E,+Er  +  ^WS{Qi-A9^)  (7) 

where  is  the  smooth  surface  emissivity,  is  the  rough  surface  and  foam  contribution  to 
the  total  emissivity,  WS  is  wind  speed  (m/s)  and  0i  is  the  incidence  angle.  The  third  term 
on  the  right  in  equation  (7)  is  the  incidence  angle  dependence  of  the  roughness  term.  The  P 
coefficient  in  this  term  contains  a  scaling  factor  between  u*  and  wind  speed  at  19.5  meters. 
The  P  coefficient  for  the  five  low  frequency  channels  of  the  SSM/I  are  in  Table  3 
(Wentz  1992). 

As  with  the  SMMR,  E^  for  the  SSM/I  frequencies  are  computed  using  the  Klein  and 

Swift  (1977)  equation  for  the  dielectric  constant  of  sea  water  and  the  Fresnel  relations 
(Smith  1988).  The  new  E^  relation  is  valid  for  incidence  angles  of51"  ±3".  The  expression 
for  Ej  at  the  five  lower  SSM/I  frequencies  is 

t^s  =  - Ys - 

where 

t=Ts-TI2,.\6  (9) 

^  =  0,--5r  (10) 

and  Ts  is  the  sea  surface  temperature.  The  values  for  the  S  coefficients  are  in  Table  3 


(Wentz  1992). 
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TABLE  3.  Coefficients  for  the  smooth  surface  emissivity  (EJ  expression  at  the 
SSM/I  channels  (Wentz  1992). 


CoefF. 

Units 

19  V 

19H 

22  V 

37  V 

37  H 

So 

K 

1.6253x10^ 

0.8220x10" 

1.6699x10" 

1.8631x10" 

0.9974x10" 

S, 

none 

-2.570x10' 

-2.805x10' 

-3.408x10' 

-5.637x10' 

-6.171x10' 

S. 

K-' 

1.729x10-^ 

1.237x10'" 

1.735x10" 

1.481x10" 

1.437x10" 

S3 

-1.177x10^ 

-0.925x10^ 

-1.036x10-" 

-0.296x10"' 

-0.707x10-" 

S4 

BUdeg 

2.162 

-1.472 

2.164 

2.123 

-1.701 

Ss 

deg-' 

0.70x10'^ 

0.21x10" 

0.75x10" 

1.17x10" 

0.55x10" 

Se 

K/deg^ 

4.5x10" 

-1.6x10'^ 

4.5x10" 

4.1x10" 

-1.9x10" 

S7 

K-'deg-* 

0.14x10-" 

-1.10x10-" 

0.02x10-" 

-0.71x10-" 

-1.27x10-" 

s/m 

4.89x10'" 

9.68x10" 

3.78x10" 

3.52x10" 

8.55x10" 

s/m 

0.623x10-" 

2.340x10" 

0.634x10" 

0.700x10" 

4.100x10" 

ni2 

s/m 

3.462x10" 

6.146x10" 

3.305x10" 

2.500x10" 

7.300x10" 

P 

s/m  deg 

-0.812x10^ 

0.806x10^ 

-0.868x10-" 

-1.193x10"' 

1.052x10"' 
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Next,  the  contribution  by  the  roughened  sea  surface  and  sea  foam  to  the  total  sea 
surface  emissivity  is  calculated.  The  brightness  temperature  due  to  the  smooth  surface  and 
the  incidence  angle  dependence  of  the  roughness  term  (third  term  in  equation  (7))  was 
subtracted  from  the  total  brightness  temperature.  Then  the  residual  brightness  temperatures 
are  fitted  by  least  squares  means  to  get  a  best  fit  equation  for  the  rough  surface  contribution 
to  the  total  emissivity.  The  resulting  expressions  are 


II 

3 

fVS  <  Imis 

(11) 

7<WS<  llm/s 

(12) 

Er  =  m2WS—  12(/n2  -m\) 

\1<WS 

(13) 

The  m  coefficients  are  in  Table  3  (Wentz  1992). 
c.  Sea  surface  emissivity  model 

Wilheit  developed  a  physically  based  sea  surface  emissivity  model.  Where  the  sea 
surface  emissivity  is  a  combination  of  the  emissivity  given  by  an  ensemble  of  flat  facets 
and  the  reduction  of  the  reflectivity  of  the  facets  due  to  a  layer  of  absorbing  nonpolarized 
foam  (Wilheit  1979). 

Lane  and  Saxton  (1952)  did  several  experiments  to  determine  the  dielectric  constant 
of  salt  water  with  varying  degrees  of  salinity  in  a  temperature  range  of  -10  to  15°C  and  a 
frequency  range  of  9  to  48  GHz.  Using  the  experimental  results  of  Lane  and  Saxton's  work 
in  conjunction  with  the  Fresnel  relations,  sea  surface  emissivity  is  computed  in  Wilheit's 


model  (Wilheit  1979). 
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Wilheit's  model  does  not  compute  the  smooth  surface  component  of  the  surface 
emissivity  separately  from  the  rough  surface  component  as  in  the  Wentz  functions.  Rather, 
the  model  treats  the  sea  surface  as  an  ensemble  of  facets  (Cox  and  Munk  1955).  Cox  and 
.  Munk  have  shown,  from  sun  glint  patterns  at  known  wind  speeds,  that  the  distribution  of 
facet  slopes  is  Gaussian.  The  distribution  that  they  used  is  expressed  as 


P(^Zz,Zy) 


(14) 


where  is  the  variance  of  the  facet  slopes.  There  is  a  factor  of  2  missing  in  the 
denominator  of  the  term  in  front  of  the  exponential  term  and  in  the  exponential  term,  which 
is  present  in  the  modem  definition  of  a  normal  distribution.  This  is,  however,  the  form 
used  by  Cox  and  Munk.  The  distribution  without  the  factor  of  2  is  narrower  than  one  with 
the  factor  of  2.  The  mean  square  slope  of  the  facets  increases  as  the  wind  speed  increases, 
so  that  cr^  be  written  as  a  function  of  wind  speed  as  follows: 


=0.003 +  5. 12  X 10-3 0.004  ^^5) 

Wind  speed  (WS)  is  at  12.5  meters  above  the  surface  (Cox  and  Munk  1955).  The  Wilheit 
model  makes  a  slight  adjustment  in  the  variance  to  account  for  a  different  wind  speed 
height  of  19.5  meters.  The  variance  in  the  model  is  given  by  the  expression  (Wilheit  1979) 

0^=0.003+4.8x10-31^5  (16) 

This  variance  of  the  sea  surface  roughness  leads  to  an  overestimate  of  the  roughness  at 
microwave  frequencies  (Hollinger  1971).  From  measurements  of  wind  speed  and 
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microwave  brightness  temperatures,  Hollinger  found  that  at  1.41,  8.36,  and  19.34  GHz  the 
variance  is  approximately  1/3,  1/2,  and  2/3  of  the  Cox  and  Munk  variance  respectively. 
Wilheit  found  that  above  35  GHz  the  full  Cox  and  Munk  variance  should  be  used.  To  take 
these  findings  into  account,  Wilheit's  surface  emissivity  model  adjusts  the  Cox  and  Munk 
variance  as  follows; 


a2(/)  =  (0.3+0.02y)aL 

f<Z5GHz 

(17) 

o^(/)  =  aL 

f>35GHz 

(18) 

where  is  the  Cox  and  Munk  variance  and  f  is  frequency  in  GHz  (Wilheit  1979). 

Having  described  the  sea  surface  as  a  Gaussian  distribution  of  facets,  the  model 
calculates  the  surface  emissivity  as  a  summation  of  each  facet's  emissivity.  Each  facet's 
emissivity  is  calculated  by  using  an  equation  derived  from  Lane  and  Saxton's  dielectric 
constant  data  and  the  Fresnel  relation;  with  each  facet  being  treated  as  a  flat  (smooth) 
surface  at  a  different  incidence  angle. 

Once  the  surface  emissivity  is  calculated,  the  model  takes  the  sea  foam  effect  into 
account.  Sea  foam  increases  the  radiance  seen  by  the  sensor.  The  increase  in  sea  surface 
emissivity  is  from  0.37  to  0.6-0.8  at  10  GHz  (Williams  1971).  During  an  experiment  over 
the  North  Atlantic,  Nordberg  et  al.  (1971)  showed  that  changes  in  brightness  temperature 
increase  linearly  with  wind  speeds  above  7  m/s.  The  amount  that  foam  affects  the 
emissivity  is  a  function  of  frequency  and  proportional  to  wind  speed  above  7  m/s  (Wilheit 
1979).  The  model  decreases  surface  reflectivity  due  to  foam  by 


E/~  a  1  -  exp  j- 


iWS-7) 


WS  >  Im/s 


(19) 
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£/=0  WS<lmls  (20) 

where  a=0.006  s/m,  and  fo=7.5  GHz.  This  function  is  the  best  fit  to  the  experimental  data 
provided  by  Webster  et  al.  (1976)  and  Shutko  (1978). 
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3.  REFBVEMENT  AND  COMPAKISON 

A  better  model  for  the  sea  surface  emissivity  is  needed  to  investigate  which 
frequencies  will  produce  the  best  set  of  brightness  temperatures  to  be  able  to  retrieve  the 
atmospheric  variables.  The  two  existing  Wentz  functions  for  surface  emissivities  are  a  best 
fit  to  available  data.  They  have  some  residual  errors  and  are  not  totally  self-consistent. 
The  actual  values  of  the  functions  may  not  be  the  same,  but  the  form  of  the  output  is 
consistent.  Therefore,  the  changes  to  Wilheit's  model  will  have  the  same  form  of  the 
Wentz  functions  but  not  necessarily  the  same  values. 

Several  changes  are  made  to  Wilheit's  sea  surface  emissivity  model.  The  first 
change  is  to  the  model's  treatment  of  multiple  reflections.  The  second  change  is  made  to 
the  sea  surface  roughness  parameter.  Finally,  the  last  change  is  in  the  treatment  of  sea 
foam. 

The  following  variables  are  held  constant  during  the  process  of  refining  Wilheit's 
sea  surface  emissivity  model.  An  arbitrary  sea  surface  temperature  of  290  K  is  chosen. 
The  SMMR  sea  surface  emissivity  function  uses  a  fixed  incidence  angle  of  49”.  Therefore, 
a  fixed  incidence  angle  of  49°  is  used.  Finally,  the  wind  speed  ranged  from  0  to  40  m/s. 
These  variables  stay  constant  throughout  the  research  unless  otherwise  specified. 

The  first  comparison  between  the  model  and  the  Wentz  SSM/I  function  for  sea 
surface  emissivity  is  the  angle  dependence  of  the  emissivity.  This  comparison  was  only 
accomplished  for  the  SSM/I  function  because  the  SMMR  function  is  only  valid  at  a  fixed 
incidence  angle.  The  SSM/I  function  is  valid  for  the  incidence  angle  range  of  48”  to  54°. 
The  incidence  angle  is  varied  between  45°  and  55°.  The  model's  incidence  angle 
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dependence  followed  the  Wentz  function  for  surface  emissivity  (see  Fig.  2).  Therefore,  the 
model  handles  the  incidence  angles  correctly. 

The  model  also  handled  sea  surface  temperatures  correctly.  In  Fig.  3,  it  is  seen  that 
emissivity  decreases  with  increasing  sea  surface  temperature.  This  is  the  same  pattern  that 
would  result  from  the  Wentz  functions  because  of  the  l/T,  relationship  in  (see  Equations 
(2)  and  (8)). 

Next,  the  wind  speed  dependence  of  the  model  is  compared  to  the  Wentz  functions. 
As  wind  speed  increases,  the  model's  emissivity  departs  more  from  the  emissivity  given  by 
the  Wentz  functions  (see  Fig.  4).  The  sharp  increase  in  emissivity  after  7  m/s  is  because  at 
this  wind  speed  the  model  turns  on  the  foam  effect.  Also,  the  Wentz  function  shows  a 
continuous  increase  in  emissivity  but  the  model  tends  to  decrease  the  emissivity  in  the 
vertical  polarization  until  the  foam  effect  starts.  The  decrease  in  the  model  emissivity 
between  0  and  7  m/s  is  because  the  foam  effect  does  not  start  until  7  m/s.  In  the  real 
atmosphere  there  is  a  smooth  transition  to  increasing  foam  that  starts  below  7  m/s.  This  is 
best  seen  in  the  37  vertical  channel  (see  Fig.  4). 

The  way  that  the  model  handles  multiple  reflections  is  adjusted  first  because  it 
affects  the  roughness.  In  the  original  model,  all  multiple  reflections  are  absorbed.  This 
yields  a  much  higher  emissivity  than  is  given  by  the  Wentz'  sea  surface  emissivity  functions 
(see  Fig.  4).  One  would  expect  that  some  of  the  background  radiation  would  go  through 
multiple  reflections  and  still  be  seen  by  the  microwave  sensor.  If  we  allow  all  the  multiple 
reflections  to  remain,  we  have  a  better  fit  to  the  emissivity  given  by  the  Wentz  functions 
(see  Fig.  5).  The  treatment  of  multiple  surface  interactions  in  this  manner  is  acceptable 
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Fig.  2.  Incidence  angle  dependence  of  the  model  and  the  SSM/I  function  for 
surface  emissivity. 
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Fig.  4.  Wind  speed  dependence  of  the  model  and  the  SSM/I  function. 
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Fig.  5.  The  treatment  of  multiple  reflections  in  the  model. 
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because  if  a  ray  is  reflected  out  of  the  field  of  view,  on  the  average,  another  ray  is  reflected 
back  into  the  field  of  view.  Also  shown  in  Fig.  5  is  the  case  where  the  multiple  reflections 
are  modeled  as  though  the  second  surface  interaction  were  at  nadir.  The  case  where  all  the 
radiation  is  reflected  back  into  the  view  path  follows  just  below  the  emissivities  given  by 
the  Wentz  functions.  Therefore,  this  case  is  used  in  the  model  for  the  rest  of  the 
calculations. 

Next,  the  roughness  parameter  of  the  model  is  compared  to  the  roughness  given  by 
the  Wentz  functions.  The  roughness  parameter  (F)  is  defined  as 


^  ~  Ry~  l-Ev 


(21) 


where  Rt  is  the  horizontal  reflectivity,  Et  is  the  horizontal  emissivity,  is  the  vertical 
reflectivity,  and  is  the  vertical  emissivity.  This  parameter  does  not  have  any  foam  effect 
because  foam  has  no  polarization  characteristics  (Wilheit  1979).  From  the  roughness  plots 
at  high  wind  speeds,  more  roughness  is  needed  at  37  GHz  and  at  frequencies  below  19.35 
GHz  less  roughness  is  needed  (see  Fig.  6). 

Roughness  in  the  model  is  adjusted  by  changing  the  variance  in  the  mean  square 
slopes  of  the  facets.  At  a  wind  speed  of  20  m/s,  the  6.6,  10.7,  18,  and  37  GHz  channels 
need  a  change  of  -10.3%,  -6.8%,  -1.3%,  and  32%  in  roughness  respectively.  This  change 
is  needed  even  after  taking  into  account  the  change  to  the  Cox  and  Munk  variance  made  by 
Hollinger  (1971).  A  line  fit  to  the  above  data  points  yield  a  decrease  in  roughness  below 
16.6  GHz  and  an  increase  above  16.6  GHz  (see  Fig.  7).  The  crossover  at  16.6  GHz 
suggests  that  Hollinger  overestimates  the  slope  variance  below  this  point  and 
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Fig.  6.  Plots  of  the  roughness  parameter  for  SMMR  frequencies. 
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Fig.  7.  A  comparison  of  Hollinger's  and  the  new  model's  Cox  and  Munk  variance 
adjustment. 
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underestimates  the  variance  above  this  point.  The  resulting  equation  is  added  to  the 
reduction  of  the  Cox  and  Munk  variance  as  found  by  Hollinger.  The  new  adjustment  to  the 
Cox  and  Munk  variance  is  given  by  the  following  expression; 

0.06852  +  0.03393  X/  (22) 

where  f  is  the  frequency  in  GHz.  This  adjustment  to  the  roughness  fits  the  data  better  and 
closer  to  the  roughness  given  by  the  Wentz  SSM/I  37  GHz  function  (see  Fig.  8).  The  new 
expression  yields  a  132%  of  the  Cox  and  Munk  variance  at  37  GHz.  Guillou  et.  al.  (1995) 
have  suggested  that  at  high  frequencies  the  sea  surface  roughness  may  be  as  high  as  2.5 
times  the  Cox  and  Munk  roughness.  At  frequencies  above  40  GHz,  the  model  uses  100%  of 
the  Cox  and  Munk  variance.  At  the  40  GHz  cutoff,  the  model  uses  143%  of  the  Cox  and 
Munk  variance.  This  cutoff  is  chosen  because  not  much  data  are  available  above  this 
frequency. 

Having  a  good  estimate  of  the  roughness,  we  can  now  make  an  adjustment  to  the 
foam  contribution.  The  emissivity  is  now  computed  without  foam  using  the  model.  This 
emissivity  is  subtracted  from  the  total  emissivity  obtained  from  the  Wentz  functions.  The 
remaining  emissivity  is  considered  to  be  from  foam.  The  residuals  should  be  zero  at  0  m/s 
wind  speed,  but  they  are  not  (see  Fig.  9).  The  model  emissivities  are  considered  to  be  more 
accurate  because  the  computed  smooth  surface  emissivities  are  based  on  the  Lane  and 
Saxton  dielectric  data  rather  than  the  Klein  and  Swift  expression  for  the  dielectric  constant. 
Therefore,  the  offset  at  Om/s  wind  speed  is  subtracted  from  all  wind  speeds  at  each 
frequency.  This  adjusts  for  zero  foam  at  0  m/s  wind  speed  (see  Fig.  10). 


Fig.  8.  Plots  of  the  roughness  parameter  for  SMMR  frequencies  after  the  variance 
adjustment. 
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The  graph  of  the  residuals  shows  two  linear  segments;  one  below  wind  speeds  of 
7  m/s  and  one  above  17  m/s.  This  result  can  be  expected  since  Wentz  developed  the  sea 
surface  emissivity  functions  in  this  manner  with  respect  to  wind  speed.  He  connected  the 
two  linear  segments  with  a  quadratic  spline.  In  the  same  manner,  a  base  function  for  the 
foam  effect  is  developed  which  is  only  a  function  of  wind  speed  (see  Fig.  11).  This 
function  is  chosen  so  that  when  it  is  subtracted  from  the  6.6  GHz  horizontal  and  vertical 
channels  the  horizontal  remainder  stays  positive  and  is  close  to  the  6.6  GHz  vertical  values. 
The  base  function  chosen  is 


O.OOUxWS 

WS<7m/s 

(23) 

0.000195  X  WS^  -0.00153  x  ^5+0.009555 

7<WS<nmls 

(24) 

0.0051  xIF5- 0.04587 

17  <WS 

(25) 

When  subtracted  from  the  total  foam  effect,  this  function  reduces  the  bend  in  the  foam 
contribution  (see  Fig.  12). 

After  the  base  function  is  subtracted  from  the  foam  contribution  at  each  frequency 
and  polarization,  the  remaining  foam  effect  is  a  function  of  frequency  and  wind  speed.  To 
determine  the  frequency  dependence,  the  remaining  foam  contribution  is  expressed  as  a 
fraction  of  the  base  function  (see  Fig.  13).  From  Fig.  13,  the  fractions  of  the  base  function 
remain  fairly  constant  below  7  m/s  and  above  17  m/s;  except  in  the  37  GHz  channel.  An 
average  fraction  is  chosen  for  each  of  these  segments  at  each  frequency  and  polarization. 
At  37  GHz,  the  curvature  to  the  fraction  is  due  to  atmospheric  interference  which  makes  the 
37  GHz  channel  difficult  to  model.  The  37  GHz  emissivities  from  the  Wentz  SMMR 


5  10  15  20  25  30  35  40 

Wind  Speed 

Horizontal  Polarization 

Fig.  13.  The  remaining  foam  contribution  expressed  as  a  fraction  of  the  base 
function  (the  straight  line  segments  are  the  new  functions  for  the  remaining  foam 
contributions). 
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and  SSM/I  functions  do  not  even  agree  (see  Fig.  14).  The  functions  are  not  totally 
self-consistent  because  they  are  a  best  fit  to  the  data.  Also,  there  does  not  exist  any  data 
that  would  suggest  that  the  37  GHz  foam  fraction  should  behave  differently  than  the  foam 
fraction  at  the  other  frequencies.  Therefore,  an  average  fraction  can  be  used  for  37  GHz 
without  much  error.  A  polynomial  is  fit  to  the  average  fraction  as  a  function  of  frequency 
for  both  horizontal  and  vertical  polarizations  (see  Fig.  13).  The  resulting  equations  are  in 
the  form  of 

aQ+a\'x.f+a2'Kf- +a^y.f  (26) 

where  f  is  frequency  (GHz).  The  a  coefficients  are  in  Table  4. 

Between  7  m/s  and  17  m/s,  the  remaining  foam  contribution  is  a  function  of 
frequency  and  wind  speed.  The  function  chosen  is  linear  in  wind  speed  and  connects  the 
constant  fractions  below  7  m/s  and  above  17  m/s  at  each  frequency  (see  Fig.  13).  The 
slopes  of  these  connecting  segments  are  different  for  each  frequency  and  polarization.  The 
slopes  and  intercepts  of  the  connecting  segments  are  calculated  for  each  frequency  and 
polarization.  Next,  polynomials  are  fit  to  the  slopes  and  intercepts  as  a  function  of 
frequency.  The  slope  and  intercept  equations  are  combined  to  form  one  equation;  one  for 
the  horizontal  polarization  and  one  for  the  vertical  polarization.  The  resulting  equations  are 
in  the  form  of 

bo  +  b\  x/-i-/>2x/  +  />3x/  +  (Co+Ci  y.f+C2'X.f)'xWS  (27) 

This  equation  is  valid  for  wind  speeds  between  7  and  17  m/s.  The  b  and  c  coefficients  are 


in  Table  4. 
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TABLE  4.  Coefficients  for  the  frequency  dependence  equations  of  the  foam  effect. 


Horizontal 

Polarization 

Vertical 

Polarization 

ao 

-0.36025 

-0.58956 

-0.7783 

-0.3146 

ai 

0.12506 

0.1478 

0.1643 

0.1325 

a2 

-0.0051709 

-0.0062256 

-0.005604 

-0.005697 

as 

0.000092046 

0.00008869 

0.0001066 

0.00008414 

l<WS<\lmls 

bo 

-0.1991 

b, 

0.109 

b. 

-0.004426 

-0.005541 

bs 

0.00009429 

0.0001223 

Co 

-0.02475 

0.03423 

Cl 

0.002654 

-0.0006405 

C2 

-0.00012682 

-0.0001519 
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The  overall  frequency  dependence  of  the  foam  effect  is  similar  to  frequency 
dependence  given  by  Wilheit's  foam  function.  His  function  is  exponential  in  frequency  and 
the  new  foam  function  is  cubic  in  frequency.  This  new  function  is  similar  to  Wilheit's 
function  until  approximately  17  GHz  at  which  point  it  begins  to  increase  (see  Fig.  15). 
This  increase  in  the  function  is  needed  to  account  for  the  larger  foam  effect  seen  at  37  GHz. 

The  overall  wind  speed  dependence  of  the  foam  effect  is  similar  to  the  base  function 
(see  Fig.  11).  Also  in  Fig.  11,  the  wind  speed  dependence  of  Wilheit's  function  is  seen. 
The  new  wind  speed  dependence  of  the  foam  effect  is  a  smooth  transition  to  increasing 
foam  as  the  wind  speed  increases.  The  new  foam  function  also  has  more  foam  at  all  wind 
speeds. 

The  resulting  sea  surface  emissivity  model  is  a  big  improvement  over  the  old 
model.  The  new  model's  emissivity  follow  closely  to  the  emissivities  given  by  the  Wentz 
functions  (see  Fig.  16).  The  emissivities  from  the  new  model  are  slightly  offset  from  the 
emissivities  given  by  the  Wentz  functions.  The  offset  is  by  design  because  the  computed 
emissivities  at  0  m/s  wind  speed  are  different.  The  emissivities  given  by  the  new  model 
follow  the  form  of  the  Wentz  functions.  The  emissivities  from  the  new  model  are  within 
2%  of  the  emissivities  from  the  Wentz  functions. 

The  fortran  source  code  for  the  new  sea  surface  emissivity  model  is  given  in 
Appendix  A.  This  code  is  used  to  produce  the  total  sea  surface  emissivity. 
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Fig  15.  Frequency  dependence  of  the  foam  effect  at  20  m/s  wind  speed 
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.  16.  Emissivities  from  the  new  model  and  the  Wentz  functions. 
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4.  CONCLUSION 

Several  changes  are  made  to  Wilheit's  sea  surface  emissivity  model.  The  first 
change  to  the  model  is  to  the  model's  treatment  of  multiple  reflections.  Multiple  reflections 
are  now  treated  as  if  the  radiation  were  reflected  back  into  the  view  path  of  the  microwave 
sensor.  This  change  lowered  the  emissivity  of  the  sea  surface;  which  is  more  representative 
of  observations  without  sea  foam.  The  second  change  is  made  to  the  sea  surface  roughness 
parameter.  An  increase  in  roughness  is  needed  at  frequencies  above  16.6  GHz  and  a 
decrease  below  16.6  GHz.  The  roughness  is  increased  to  132%  of  the  Cox  and  Munk 
roughness  at  37  GHz  and  30%  of  the  roughness  at  6.6  GHz.  The  last  change  to  the  model 
is  in  the  treatment  of  sea  foam.  The  foam  effect  is  now  a  smooth  transition  to  increasing 
foam  as  the  wind  speed  increases;  instead  of  being  switched  on  at  7  m/s.  These  changes 
yield  an  improved  sea  surface  model  whose  emissivities  fall  within  2%  of  the  emissivities 
from  the  Wentz  functions. 

While  the  Wentz  sea  surface  emissivity  functions  are  adequate  for  the  SMMR  and 
the  SSM/I  frequencies  and  incidence  angles,  there  remains  a  need  for  a  more  general  model 
for  all  frequencies  and  incidence  angles.  The  new  model  does  not  have  these  limitations. 
Therefore,  it  can  be  used  to  simulate  brightness  temperatures  for  other  frequencies  and 


incidence  angles. 


38 


REFERENCES 

Chang,  A.  T.  C.,  and  T.  T.  Wilheit,  1979;  Remote  sensing  of  atmospheric  water  vapor, 

liquid  water,  and  wind  speed  at  the  ocean  surface  by  passive  microwave  techniques 
from  the  Nimbus  5  satellite.  Radio  Sci.,  14,  793-802. 

Cox,  C.,  and  W.  Munk,  1955:  Some  problems  in  optical  oceanography.  J.  Marine  Res.,  14, 
63-78. 


Griffiths,  D.  J.,  1989:  Introduction  to  Electromagnetics  (2nd Ed).  Prentice  Hall,  532  pp. 

Guillou,  C.,  S.  J.  English,  C.  Prigent,  and  D.  C.  Jones,  1995:  Passive  microwave  airborne 
measurements  of  the  sea  surface  response  at  89  and  157  GHz.  Submitted  to 
J.  Geophys.  Res.  Personal  collection  of,  D.  J.  Kohn. 

Hollinger,  J.  P.,  1971:  Passive  microwave  measurements  of  sea  surface  roughness.  IEEE 
Trans.  Geosci.  and  Remote  Sensing,  GE-9,  165-169. 

Klein,  L.  A.,  and  C.  T.  Swift,  1977:  An  improved  model  for  the  dielectric  constant  of  sea 
water  at  microwave  frequencies.  IEEE  Trans.  Antennas  Propagat.,  AP-25, 

104-111. 

Lane,  J.  A.,  and  J.  A.  Saxton,  1952:  Dielectric  dispersion  in  pure  polar  liquids  at  very  high 
radio-frequencies  I.  Measurements  on  water,  methyl  and  ethyl  alcohols.  Proc.  Roy. 
Soc.,  London  A,  214,  400-408. 

Nordberg,  W.,  J.  Conaway,  D.  B.  Ross,  and  T.  Wilheit,  1971;  Measurements  of  microwave 
emission  from  a  foam-covered,  wind-driven  sea.  J.  Atmos.  Sci.,  28, 429-435. 

Ross,  D.  B.,  V.  Cardone,  1974:  Observations  of  oceanic  whitecaps  and  their  relation  to 
remote  measurements  of  surface  wind  speed.  J.  Geophys.  Res.,  79,  444-452. 

Shutko,  A.,  1978;  Report  on  Soviet  progress  in  microwave  radiometry  of  the  ocean's 
surface.  lUCRM  Colloquium  on  Passive  Radiometry  of  the  Ocean's  Surface. 

Patricia  Bay,  B.  C.,  Canada. 

Smith,  P.  M.,  1988;  The  emissivity  of  sea  foam  at  19  and  37  GHz.  IEEE  Trans.  Geosci. 
and  Remote  Sensing,  26,  541-547. 

Wentz,  F.  J.,  1983:  A  model  function  for  ocean  microwave  brightness  temperatures.  J. 
Geophys.  Res.,  88,  1892-1908. 


39 


_ ,  L.  A.  Mattox,  and  S.  Peteherych,  1986:  New  algorithms  for  microwave 

measurements  of  ocean  winds  applications  to  SEASAT  and  the  special  sensor 
microwave  imager.  J.  Geophys.  Res.,  91,  2289-2307. 

_ ,  1992:  Measurement  of  oceanic  wind  vector  using  satellite  microwave  radiometers. 

IEEE  Tram.  Geosci.  and  Remote  Sensing,  30,  960-972. 

Wilheit,  T.  T.,  1979:  A  model  for  the  microwave  emissivity  of  the  ocean's  surface  as  a 
function  of  wind  speed.  TREE  Tram.  Geosci.  Electron.,  G1E-\1, 244-249. 

_ ,  and  A.  T.  C.  Chang,  1980:  An  algorithm  of  ocean  surface  and  atmospheric 

parameters  from  the  observations  of  the  scaiming  multichannel  microwave 
radiometer.  Radio  Sci.,  15,  525-544. 

Williams,  G.  F.,  1971:  Microwave  emissivity  measurements  of  bubbles  and  foam.  IEEE 
Tram.  Geosci.  Electron.,  GE-9,  221-224. 


n  o  n  o  o 


40 


APPENDIX  A 

SUBROUTINE  ROUGH(F, ANGLE, TEMP, PPT.WS  ,VERT,HORIZ) 

c 

C  ♦  TfflS  SUBROUTINE  COMPUTES  VERTICAL  AND  HORIZONTAL  * 

C  *  EMISSrVITIES  FROM  COX  AND  MONK  AT  20M  WHERE  * 

C  *  * 

C  *  F-FREQUENCY  IN  GHz  * 

C  *  ANGLE-INCIDENCE  ANGLE  IN  RADS  * 

C  *  TEMP-TEMP  INK  * 

C  *  PPT-SALINITY  IN  PARTS  PER  THOU.  * 

C  *  WS-WIND  SPEED  IN  M/S  * 

C  *  VERT- VERTICAL  POLARIZATION  * 

C  *  HORIZ-HORIZONTAL  . .  * 

C 

C  Variance  of  sea  surface  roughness  given  by  Cox  and  Munk.  And  changes  to  the 
C  variance  as  given  by  the  new  model. 

C 

VAR=0.003+4.80*WS/1000. 

IF  (F  .LE.  40.0)  VAR=VAR*(0.06852-K).03393*F) 

T=TEMP 
VERT=0. 

HORIZ=0. 

YNORM=0. 

SIG=SQRT(VAR) 

STEP=SIG/10. 

SA=SIN(ANGLE) 

CA=COS(ANGLE) 

Calculation  of  the  relaxation  time  (TAU)  of  the  dielectric  constant 
of  sea  water  (Lane  and  Saxton) 
and  International  Critical  Tables  conductivity  data 

S=PPT/58.4 

XLAM=5 1  .-4.83  5*S+1 .728*(T-273  .)-0.2037*S*(T-273 .) 

IF(S.GT.  1.)  XLAM=XLAM-0.35/S 
ES=190.-81*S+38*S*S-(03.75-2.*S+S*S)*T/10. 

T1=0.00199*EXP(2140/T)/T 
T2=0.00243*EXP(2060fr)/T-Tl 
T3=0.00324*EXP(1 968/T)/T-T  1 
TAU=T1+(4.*T2-T3)*S-K2.*T3-4.*T2)*S*S 
C 


non  noon  non  non  ooooo 
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C  TAU  in  nanoseconds,  F  in  GHz,  T  in  deg  K,  S  is  salinity 
C  Dielectric  constant  Calculation 

C  RPE,  XIPE  are  the  real  and  imaginary  parts  of  the  dielectric  constant  respectively 
C 

X=6.283*F*TAU 
EZ=4.9 

RPE=EZ+(ES-EZ)/(1  .+X*X) 

XIPE=(ES-EZ)*X/(1  .+X*X)+1 .8*XLAM*S/F 

Taking  the  square  root  of  the  complex  dielectric  constant 
to  get  the  complex  index  of  refraction 

RPN,  XEPN  are  the  real  and  imaginary  parts  of  the  index  of  refraction  respectively 

XMODE=SQRT(RPE**2+XIPE**2) 

XMODN=SQRT(XMODE) 

ARGN=0.5*ATAISr(XIPE/RPE) 

RPN=XMODN*COS(ARGN) 

XIPN=XMODN*SIN(ARGN) 

Calculates  the  emissivity  contribution  of  each  facet  and  sums  up  the  emissivities 

DO  1001=1,61 
DO  100KK=1,31 
K=KK-1 

ZX=STEP*(J-31) 

ZY=STEP*K 


The  distribution  used  the  one  used  by  Cox  and  Munk 

WEIGHT=EXP(-((J-3 1)**2+K*K)/100.) 
IF(KK.EQ.l)  WEIGHT=WEIGHT/2. 


Facet  hidden  from  view 

RDN  is  the  dot  product  between  the  incoming  ray  and  facet  normal 

RDN=CA+ZX*SA 
IF  (RDN  XT.  0.)  GO  TO  100 


STP  is  the  sine  squared  of  above  angle  calculated  by  1 -cosine*  *2 

STP=(SA*SA+ZX*ZX+ZY*ZY-(ZX*SA)**2-2.*ZX*SA*CA)/(1+ZX*ZX+ZY*ZY) 

SI=SQRT(STP) 

IF  (SI  .LT.  0.001)  SI=0.001 
EM=ZY*ZY/(1  .+ZX*ZX+ZY*ZY) 


ooo  oiro  o  ooo  noooo 
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EM=EM/(SI*SI) 

IF  (SI  XT.  O.OOIOl)  EM=  5 
CI=SQRT(1.-SI*SI) 
WEIGHT=WEIGHT*CI 
YNORM=YNORM+WEIGHT 


Projection  of  facet  in  view  direction 

Squaring  the  complex  index  of  refraction  and  finds  the  complex  conjugate 
of  the  squared  index  of  refraction 


AN2=RPN**2+XIPN**2 

SRR=RPN*SI/AN2 

SRI=-XIPN*SI/AN2 

XMOEK(1-SRR**2+SRI**2)**2-K2.*SRR*SRI)**2)**0.25 

ARG=-0.5*ATAN(2.*SRR*SRI/(1.-SRR**2+SRI**2)) 

CRR=XMOD*COS(ARG) 

CRI=XMOD*SIN(ARG) 


Calculates  the  reflected  emissivity  with  the  Fresnel  equations 


RH=((SI*CRR-CI*SRR)**2+(SI*CRI-CI*SRI)**2)/((SI*CRR+CI*SRR)**2+ 

&(SI*CRI+SRI*CI)**2) 

RV=RH*((CI*CRR-SI*SRR)**2-KCI*CRI-SI*SRI)**2)/((CI*CRR+SI*SRR)**2+ 

&(CI*CRI+SI*SRI)**2) 

HORIZ=HORIZ+WEIGHT*((l  .-EM)*(1  .-RH)+EM*(1  .-RV)) 
VERT=VERT+WEIGHT*((1  .-EM)*(1  .-RV)+EM*(1  .-RH)) 

)0  CONTINUE 

HORIZ=HORIZATSfORM 

VERT=VERTAT^ORM 


Adds  foam  effect 

IF  (WS  .LE.  7.0)  THEN 
C 

FOAM=0.0012*WS 

FOAMV=FOAM+FOAM*(-0.7783+0. 1 643  ♦F-0.005604*F*F-K).0001 066*F*F*F) 
FOAMH=FOAM+FOAM*(-O.36025-K).  12506*F-0.005 1709*F*F+ 
&0.000092046*F*F*F) 

C 

ELSEBF  ((WS  GT.  7.0)  AND.  (WS  .LT.  17.0))  THEN 
C 
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FOAM=0.000195*WS*WS-0.00153*WS+0.009555 
F0AMV=F0AM+F0AM*(-1 . 103+0. 1866*F-0.005541  *F*F+ 
&0.0001223*F*F*F+(0.03423-0.0006405*F-0.000141519*F*F)*WS) 
FOAMH=FOAM+FOAM*(-0. 1991+0. 109*F-0.004426*F*F+ 
&0.00009429*F*F*F+(-0.02475+0.002654*F-0.00012682»F*F)*WS) 

C 

ELSE 

C 

FOAM=0.005 1  *WS-0.04587 

FOAMV=FOAM+FOAM*(-0.3146+0.1325*F-0.005697*F*F+0.00008418*F*F*F) 

FOAMH=FOAM+FOAM*(-0.58956+0.1478*F-0.0062256*F*F+0.00008869*F*F*F) 

C 

ENDIF 

C 

HORIZ=FOAMH+(l.-FOAMH)*HORIZ 
VERT  =FOAMV+(l  .-FOAMV)*VERT 
RETURN 
END 
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